Divergent Signal-to-Noise Ratio and Stochastic Resonance in Monostable Systems 
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We present a class of systems for which the signal-to-noise ratio always increases when increasing 
the noise and diverges at infinite noise level. This new phenomenon is a direct consequence of 
the existence of a scaling law for the signal-to-noise ratio and implies the appearance of stochastic 
resonance in some monostable systems. We outline applications of our results to a wide variety of 
systems pertaining to different scientific areas. Two particular examples are discussed in detail. 

PACS numbers: 05.40. -fj 

Stochastic Resonance (SR) is a phenomenon wherein the response of a system to a driven periodic signal 

is enhanced at an optimized non-zero noise level. Although, increasing the noise level in order to enable us to more 
— ' easily detect a signal was considered counterintuitive, this constitutes one of the most surprising results of the SR. 
^ , It seems obvious, however, that the signal-to-noise ratio (SNR) must go to zero as noise is increased indefinitely. 
' Contrarily, in this letter we present a class of systems in which the SNR always increases when the noise is increased 
and diverges at infinite noise level, instead of exhibiting a maximum at a particular value of the noise. This result 
implies the presence of SR in monostable systems for which a maximum in the SNR at non-zero noise level has never 
T-H , been observed before. These findings open up new possibilities concerning the application of SR to a great variety of 
K>" ' physical, chemical and biological systems. To be explicit we have applied our results to two particular cases; namely, 

a ferromagnetic particle and a standard model of neural excitable medium. 
^— ^ ' The class of systems we will discuss are described by only one relevant degree of freedom whose dynamics is governed 
. by the following Langevin equation 

= (1) 

■ where h{t) = fc(l -f asin(a;o^)), with k and a (< 1) constants, n is an integer number and S^it) is Gaussian white 

^ ' noise with zero mean and second moment {£,[t)£,{t -\- t)) = DS{t), defining the noise level D. The system can be 

I characterized by the quantity v which is a function of the variable x(t). This quantity is sometimes referred to as the 

' O response to the oscillating force. The effect of this force may be analyzed by the power spectrum 
Ch , 

I p27r/ujQ poo 

O' P{^)= dt {v{t)v{t + t)) e-"^-" dT . (2) 
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To this purpose we will assume that it consists of a delta function centered at the frequency ojq plus a function Q{ijj) 
which is smooth in the neighborhood of loq and is given by P{uj) = Q{i^) + S{ujq)5{ijJ — loq). Then, the SNR is defined 
by SNR = S'(a;o)/Q(wo) and consequently, has dimensions of inverse of time. 

The existence of a characteristic time r in our system will enable us to propose the form of the SNR through the 
simple scaling law 

SNR = /(a,woT)r-i , (3) 

where — £)"/(!+») and f{a,uJoT) is a dimensionless function, provided that v{x) does not introduce 
another characteristic time. We will suppose that for a given value of r the limit of SNR when luq goes to zero exists. 
As such, the following expression for small driving frequencies holds 

SNR« /(a,0)T-\ (4) 

Let us now discuss the main characteristics of our model upon varying the exponent n. If n = 0, one finds the exact 
result SNR = f(a,uJok~^)k, which does not depend on the noise level. Even more interesting is the behavior obtained 
for the case n > 0. The scaling of the SNR indicates that it increases when increasing the noise level, achieving the 
behavior SNR cx D as n goes to infinity. A particular and common situation illustrating this case corresponds to the 
potential Vi{x) = jh{t)x'^ (Fig- obtained when n = 1, for which SNR increases as V^. In Fig. 0b we have 
depicted the SNR corresponding to Vi{x). Here the magnitude giving the response of the system has been taken 
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v{x) = x^. Our result is obtained from numerical simulations by integrating the corresponding Langevin equation by 
means of a standard second-order Runge-Kutta method for stochastic differential equations [ p^|JT^ . In order to verify 
the scaling law proposed through eq. (jj) we have fit the values of the SNR to a power law in the range of D from 1 
to 1000 for the potential Vi{x). SNR is given by aD'' with b = 0.498 ± 0.004 and a = 0.76 ± 0.02, which is in good 
agreement with the theoretical value 6 = ^. The power spectrum corresponding to Vi{x) for two values of the noise 
levels is shown in Fig. |l|c. Both the signal and noise background increase when increasing D, but the signal increases 
faster. 

In spite of its simplicity, our model encompasses a great variety of common situations. Around an equilibrium 
state most systems may be approximated by a parabolic potential. Thus, for n = our model describes a system 
around an equilibrium state in a force field whose intensity varies periodically in time. A physical realization of 
such a system could be a dipole under an oscillating field. For n = 0, however, the SNR is independent of D. 
In order to understand the behavior of the SNR, as D is increased we must take into account the corrections to 
the parabolic approximation. Commonly, these corrections are given by a term proportional to x*. Two particular 
realizations of this situation have been analyzed numerically resulting in an increase of the SNR (Fig. |b) with the 
noise level. For V2{x) = h{t){^x'^ + jx*) (Fig. |^a), which basically corresponds to a potential that around the 
minimum grows faster than a parabolic one, we expect that the SNR is an increasing function of the noise, since for 
low noise level the potential behaves as ^h{t)x'^, whereas for high noise level as jh{t)x'^. A slightly different potential 
is V3(x) — ^x^ + jh{t)x'^ (Fig- which only differs from the previous one in the behavior at low noise level. Since 
under this circumstance this potential reduces to ^x^ and consequently it is not modulated by h{t), SNR goes to zero 
for low D. For n = 1, our model describes the dynamics of a system at the critical point of both the pitchfork and 
Hopf bifurcations | |T^ , occurring in many systems including, to mention just a few [ p"5|jT^ ], chemical reactions, models 
of populations, convection in liquids, lasers and instabilities in semiconductors. 

The divergence of the SNR is due to the fact that the potential is unbounded. It is obvious that for a bounded 
potential the noise can completely destroy the response of the system. An important consequence follows from the 
previous results: if the SNR grows for low noise level (when the potential around the minimum can be approximated 
by a potential like V2 or V^) and the potential is bounded (which implies that SNR goes to zero for large noise) then the 
SNR exhibits a maximum, thus indicating the appearance of SR. To illustrate this point we report results of numerical 
simulations for the dynamics of a ferromagnetic particle p^ , p^ under an external magnetic field and with energy of 
anisotropy h{t)sin'^9, with 6 being the angle between the magnetic moment and the axis of easy magnetization. The 
external magnetic field is then applied in the direction of the easy axis of magnetization and its intensity is as high 
as the system becomes monostable. The dynamics of the magnetization may be described by the Langevin equation 

d9 

— ^ -sine - h{t) sin^ e cose + ^(t) , (5) 
dt 

where the first term on the right hand side accounts for the interaction with an external magnetic field, the second for 
the anisotropic effects and the third is a noise source due to a random field or to thermal fluctuations. The parameter 
h{t) is assumed to be oscillatory, the reason being, for example, that the presence of oscillations of the pressure of the 
medium surrounding the particle. The response of the system is now given by cos 0, i.e. by the magnetization. For 
low noise, the potential (Fig. |^a) and the magnetization reduces to ^e^ + jh{t)6^ and 1 — ^6*^, respectively. Therefore 
our previous results apply to this case. The corresponding SNR is shown in Fig. ^ and exhibits a maximum at a 
finite noise level. This result clearly shows the existence of an optimal noise level for which the system is more sensible 
to periodic changes of the environment. In Fig. ^ we have represented the power spectrum for some values of the 
noise level. 

The next example to be considered corresponds to a standard model of a neural excitable medium. This model 
characterizes the activity generated in a slab of neural tissue comprising a very large number of closely packed and 
coupled nerve cells [p^|-^. We will consider the case of all-to-all connectivity in which spatial dependence may be 
ignored The model of neural excitable medium is given by the following equation |]l9| 

C^ = -R-'U + ^{U) + P , (6) 

describing the dynamics of the spatial average of the transmembrane potential U. Here C is the membrane capacitance, 
R the membrane resistance and P is an external current applied to the net. The nonlinear term <j>{U) is proportional 
to the gain of the neuron and accounts for its mean firing rate. Its form is usually taken to be a sigmoidal, e.g. 

0(C/)=e(l + e-''(^"^))-i , (7) 
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where v is a constant, fixing the sensitivity to excitation of the population, 9 is the threshold mean voltage and e 
a parameter depending on the structure of the net and the characteristics of the neuron. We will consider that the 
external current applied to the net fluctuates and that it may be approximated by a Gaussian white noise {{P{t)) = 
and {P{t)P{t + r)) = DS{t)) The potential function (Fig. ||a) corresponding to the variable U is given by 

V{U)^^i^lR-'u'-lln{l + e-^^^-<^^)^ . (8) 

The membrane resistance R can be modified through small changes in the permeability of a suitable ion. Let us 
consider that this modification is periodic in time. There exists a range of parameters for which the neurons of the 
net are not excited, i.e. U = 0, consequently this state is not affected by small variations of R. The situation changes 
drastically with the addition of noise. For small values of the noise, when the term (j)(U) is irrelevant, the dynamics of 
this model may be described by eq. (|^), with n — 0. In this case, the SNR for v — U"^ is independent of D, but in this 
system the interesting variable is U instead of [/^. For small noise values, the potential is symmetric in U, thus giving 
a zero SNR for the variable U. When increasing the noise level the symmetry is lost, since positive fluctuations of U 
may cause firing of the neurons. Variations of R then modulates the amplitude of the fluctuations of U, giving rise to 
periodicity in the firing rate of the neurons, which implies a non- vanishing SNR. Further increasing of the noise level 
then leads to the restoration of the symmetry. Consequently, this model exhibits SR as has been shown in Fig. 
It is interesting to remark that, in contrast with previous results concerning the appearance of SR for a single neuron 
1^,^, our results refer to the prediction of the phenomenon for an ensemble of neurons. 

Our analysis has revealed the presence of SR in a wide variety of situations, embracing different scientific areas, 
which have not been considered up to now. A methodological aspect to be emphasized, and that could be used 
in subsequent studies, is that arguments as simple as scaling laws or considerations about the symmetry may help 
us in predicting the enhancement of signals via SR. Our work, then, offers new perspectives on what concerns the 
consideration of SR as a general phenomenon that might apply to diverse systems. 
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FIG. 1. (a) Potentials Vi{x), V2{x) and V3{x), for the maximum value of h{t) (solid line) and for the minimum value 
(dashed line). Here k = 1, a = 0.5 and uJo/2-w = 0.1. (b) Behavior of the SNR for the three potentials presented previously, 
(c) Power spectrum corresponding to Vi{x) for D = 0.1 and D = 1000. 

FIG. 2. (a) Potential energy of the ferromagnetic particle (eq. (5)) as function of 6 for the maximum value of h{t) (solid 
line) and for the minimum (dashed line). The parameters are taken k = 0.3, a = 3/2 and ujo/in = 0.1. (b) SNR for the 
previous values of the parameters obtained through computer simulations, (c) Power spectrum for D = 0.05, 0.15, 0.35 and 1. 

FIG. 3. (a) Representation of the potential function V{U) for the maximum value of = h{t) (solid lino) and for the 
minimum value (dashed line). The parameters used here are C = 1, 9 = 2, e = 1, v = 10 and the resistance = h{t), with 
k = 2, a = 0.5 and uJo/2-w = 0.01 (b) SNR for the previous parameters. 
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